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ABSTRACT 

Rolling-element bearing forces vary nonlinearly with bearing deflection. Thus an ac- 
curate rotordynamic analysis requires that bearing forces corresponding to the actual 
bearing deflection be utilized. For this work bearing forces were calculated by COBRA- 
AHS, a recently developed rolling-element bearing analysis code. Bearing stiffness was 
found to be a strong function of bearing deflection, with higher deflection producing 
markedly higher stiffness. Curves fitted to the bearing data for a range of speeds and 
loads were supplied to a flexible rotor unbalance response analysis. The rotordynamic 
analysis showed that vibration response varied nonlinearly with the amount of rotor im- 
balance. Moreover, the increase in stiffness as critical speeds were approached caused a 
large increase in rotor and bearing vibration amplitude over part of the speed range com- 
pared to the case of constant bearing stiffness. Regions of bistable operation were possi- 
ble, in which the amplitude at a given speed was much larger during rotor acceleration 
than during deceleration. A moderate amount of damping will eliminate the bistable 
region, but this damping is not inherent in ball bearings. 


INTRODUCTION 

Rotordynamic response of all but very flexible rotors depends strongly on bearing 
properties. When bearings are nonlinear, accurate rotor response calculations require use 
of variable bearing properties that reflect the precise conditions encountered, rather than 
average properties. While fluid film bearings are often reasonably linear for small deflec- 
tions (although there is usually a strong speed dependence) rolling-element bearings have 
a much less linear force-displacement relationship. 

Analyses of load-stress relationships for rolling-element bearings were published by 
Lundberg and Palmgren in 1947 [1]. These relationships were not easily usable until the 
analytical additions of Jones in 1960 [2] led to the marketing of his computer code. 
Jones’ code is still the most widely used rolling-element bearing analysis tool despite its 
age. Poplawski et al. [3] recently developed the code COBRA-AHS (Computer 
Optimized Ball and Roller Bearing Analysis-Advanced High Speed) under a NASA 
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Small Business Innovative Research (SBIR) contract. It is intended to be a major 
improvement over earlier rolling-element bearing analytical tools. Up-to-date stress/life 
data, accurate stress calculations considering bearing installation press fits, and an 
interactive front end are some of its features. It accounts for 5 degrees of freedom in the 
bearing, and can calculate load-displacement data for high speed radial and angular 
contact ball bearings and also for cylindrical and tapered roller bearings. 

In a previous paper by the authors [4], COBRA-AHS was used in conjunction with a 
transient rotordynamics code to calculate response to a suddenly applied large imbalance 
simulating a blade loss. The results obtained with nonlinear bearings could not be dupli- 
cated with any average bearing stiffness value. In the present paper, nonlinear bearing 
data from COBRA-AHS is used with a steady-state unbalance response code to calculate 
rotor behavior over a range of speeds, and at imbalance levels representative of normally- 
operating aerospace turbomachinery. 

Steady-state rotordynamic response codes have been in use for a number of years. One 
of the early and still viable codes was formulated by Lund [5], It uses the transfer matrix 
method to calculate unbalance response of a flexible rotor in asymmetric bearings. This 
code was simplified for the case of symmetric bearings by Kirk [6]. Kirk’s code was fur- 
ther modified for the present work to enable the use of nonlinear bearings by iterating on 
the rotor amplitude. 


ANALYTICAL SYSTEM AND PROCEDURE 


Figure 1 is a drawing of the shaft system. It depicts a fairly stiff shaft with concen- 
trated masses (which may represent compressor or turbine wheels) at stations 5 and 7. In 


Station 

12345 6 789 10 11 



0 10 20 30 


total, 11 stations and 10 elements were 
used in the model. Radial (deep groove) 
ball bearings of 25 mm bore diameter are 
at stations 2 and 10. The shaft material 
is steel; total mass of the shaft system is 
4.6 kg. 


Length, cm 

Figure 1. Sketch of rotor 



Figure 2. Critical speed map 


Figure 2 shows the first two system 
critical speeds as a function of bearing 
stiffness. For stiffness up to about 
200 MN/m, critical speeds rise rapidly 
with stiffness, indicating significant bear- 
ing participation in the rotor motion. For 
higher stiffness values, critical speeds do 
not increase as rapidly, indicating that 
increasing amounts of motion are due to 
shaft bending rather than bearing deflec- 
tion. Mode shapes calculated at critical 
speeds, shown in figure 3 for two values 
of bearing stiffness, confirm this. For a 
bearing stiffness of 18 MN/m, most of 
the elastic deflection takes place in the 
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bearings. At the higher stiffness of 180 MN/m, on the other hand, nearly all of the elastic 
deflection is in the shaft at the first critical speed, while at the second critical speed there 
is still large bearing deflection as well as noticeable shaft bending. 



(a) Bearing stiffness = 18 MN/m 



Figure 3. Critical speed mode shapes 

COBRA-AHS was used to generate load versus deflection data for the bearings at 
speeds from 10 000 to 80 000 rpm and loads from 44 to 8800 N. These data were con- 
verted to effective stiffness values by dividing the load by the deflection. A power series 
curve of the form K = a + b e was then fitted to the stiffness data for each speed, where 

K is bearing stiffness, e is bearing 
deflection, and a, b, and c are coeffi- 
cients. Cubic polynomials were fitted to 
the coefficients of the first curve fit to 
account for the coefficients’ variation 
with speed. The coefficients resulting 
from this second curve fit were then 
supplied to the rotor response code. 
Figure 4 shows stiffness, that is, load 
divided by deflection, for several 
speeds. Both the bearing data and the 
fitted curves are shown; the curves fit 
the data very well, verifying that the 
mathematical representation of the code 
output is adequate. Note that, for all 
speeds, there is an order of magnitude change over the range of the figure. This illustrates 
the significant nonlinearity of the ball bearing used, and is representative of rolling- 
element bearings. Stiffness decreases as speed increases. This occurs because centrifugal 
loads cause a deflection within the ball to outer race contact which adds to the deflection 
due to the applied load; more deflection for the same applied load means lower stiffness. 
The bearing modeled was not preloaded, but had nearly zero clearance at zero speed; the 
outward movement of the balls causes the bearing to develop a clearance as speed 



Figure 4. Ball bearing stiffness 
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increases. For the present study, this clearance has negligible effect. However, for a 
machine subject to external vibration any clearance is undesirable. Thus ball bearings are 
commonly given an axial preload to prevent any clearance. 

Balancing criteria for rigid rotors are given by international standard ISO 1940 for 
various classes of rotating machinery (see, e.g., [7]). These criteria are also commonly 
used for moderately flexible rotors. The recommended balancing class for aircraft gas 
turbine assemblies is G6.3. Permissible imbalance is inversely proportional to operating 
speed. The criteria are conservative, recognizing that considerable balance degradation 
may occur in service. Imbalance levels studied herein are approximately 4 to 32 times 
those allowed for the example rotor at 40 000 rpm, resulting in imbalances of 3 to 
23 gram centimeters. The imbalances were applied to the rotor mass at station 7. 
Response at stations 7 and 10 (the imbalance and right bearing locations) was examined 
for speeds of 20 000 to 80 000 rpm. Although larger than allowed by the standard for a 
well-balanced rotor, these imbalances are considerably less than those produced by a 

blade loss; Fleming and Poplawski 
[4] used imbalances of 36 to 
576 gram centimeters in their study 
of transient vibration. 

RESULTS 

Synchronous vibration ampli- 
tude at bearing station 10 is shown 
in figure 5 for imbalances of 6, 12, 
and 23 gram centimeters. The most 
notable feature of this figure is the 
portrayal of bistable operation: at 
some speeds, amplitude can assume 
one of two values, depending on 
the rotor history. Arrows on the 
curves in figures 5-8 and 10-13 
show whether speed was increasing 
or decreasing when the infonnation 
was generated. For discussion pur- 
poses we will focus on the data for 
12 g cm imbalance (dotted lines in 
the figures). As speed rises above 
30 000 rpm, vibration amplitude 
increases as is usual when running 
up to a critical speed. This results in 
an increase in bearing load, shown 
in figure 6. Since ball bearing stiff- 
ness increases with amplitude, there 
is also an increase in bearing stiff- 
ness as speed increases (fig. 7). 
This in turn raises the critical speed 
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(fig. 2), so operation remains subcritical longer than would occur with constant stiffness 
bearings. Rotor bending increases with stiffness (and thus with speed), as can be seen by 
comparing figure 3(a) with figure 3(b). Finally, at 48 000 rpm the critical speed is passed 
even at the large bearing stiffness. Supercritical operation begins suddenly as the ampli- 
tude jumps down to a low value. For this case, operation is now also above the second 
critical speed, as can be discerned from figure 2 and the bearing stiffness plotted in 
figure 7. Further increases in speed produce an increase in amplitude as the third critical 
speed (above 200 000 rpm for all bearing stiffnesses investigated) is approached. 



When speed is decreased from 
the maximum, the supercritical 
mode shape results in a low bear- 
ing stiffness down to 37 000 rpm 
where the amplitude jumps up to 
rejoin the curve for increasing 
speed. 

Vibration amplitude for station 
7 (the location of the imbalance) is 
shown in figure 8. The appearance 
of this figure is similar to figure 5, 
but the amplitudes are higher as 
the critical speeds are approached 
because of rotor bending. Mode 
shape plots in figure 9 illustrate 
this. In this figure, mode shapes at 
43 000 rpm are shown for both in- 
creasing and decreasing speed. 
Because the rotor mode is a com- 
bination of the first and second 
critical speed mode shapes, the 
rotor shape is not strictly planar. 
This is shown in the figure by plot- 
ting two curves for each mode; the 
first is the amplitude in the direc- 
tion of maximum displacement 

(labeled direct), and the second is 
at right angles to this (labeled 

quadrature). For increasing speed, 
with operation below the first criti- 
cal speed and high bearing stiff- 
ness, there is considerable rotor 
bending. The quadrature compo- 
nent is quite small, indicating a nearly planar mode shape as the first critical speed mode 
dominates. For decreasing speed, the mode is nearly a straight line with amplitude much 

lower than for increasing speed. In this case, however, the quadrature component is 

significant, indicating a combination of modes. 



Speed, rpm 

Figure 8. Station 7 amplitude 
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In each of figures 5-8, the data 
for 6 g cm imbalance (solid lines) 
indicate not one, but two regions 
of bistable operation. For this 
imbalance case, the first and sec- 
ond critical speeds are traversed 
one at a time, rather than simulta- 
neously as occurs with higher 
levels of imbalance. 

Bistable operation has been 
observed both analytically and 
experimentally on rotors supported 
by squeeze film dampers (e.g., 
[8]). The reasons it occurs are the 
same as discussed above, since 
squeeze film damper stiffness 
increases with vibration amplitude similar to the behavior of a ball bearing. Indeed, 
bistable operation is a well-recognized phenomenon of nonlinear systems [9], However, 
to the authors’ knowledge, bistable operation for rotors supported by nominally 
undamped ball bearings has not been reported. This is perhaps because most high-speed 
rotors, particularly those that must traverse critical speeds, are supported on damped 
bearings. 

Comparison with constant-stiffness bearings. - Figure 10 shows vibration amplitude at 
bearing station 10 for an imbalance of 12 g cm. Two fixed values of stiffness (35 and 
105 MN/m, which nearly bracket the stiffnesses shown in figure 7) were studied; the 

variable stiffness case of figure 5 is 
also shown for comparison. The 
curves for constant stiffness show 
classical rotor behavior, with a 
degree of symmetry around the 
peak amplitude. At the lower 
(35 MN/m) stiffness, two critical 
speeds are traversed. For the 
higher 105 MN/m stiffness, only 
one critical speed appears below 
80 000 rpm. This is consistent with 
the critical speed map of figure 2. 
The variable stiffness case 

20000 30000 40000 50000 60000 70000 80000 .... , . , , , , 

Speed rpm exhibits higher amplitude them the 

Figure 10. Station 10 amplitude for fixed and 105 MN/m case over most of the 

variable bearing stiffness; imbalance = 12 g cm speed range during run up. During 

run down, however, critical speed 
passage does not occur until the speed has dropped to 37 000 rpm, and peak amplitude 
is much lower than during run up because supercritical operation is maintained until 
jump up. 




Figure 9. Rotor mode shapes. 
43 000 rpm, 12 g cm imbalance 
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Further insight into rotor 
behavior may be obtained by 
examination of rotor phase angles, 
defined as the angle by which 
the imbalance force leads the 
displacement (this definition is 
equivalent to the more common 
definition for the present case of a 
single unbalance mass, and makes 
the results easier to understand). 
Figure 1 1 plots the phase angle at 
bearing station 10 for the same 
conditions as figure 10. For all bearing stiffness cases, phase angle is initially small but 
increases as the first critical speed is approached. Total phase change for the fixed stiff- 
ness cases is close to the 180 degrees predicted in classical rotor studies. When variable 
stiffness from COBRA-AHS is used, phase change is considerably less, but this may be 
explained by recalling that the rotor remains subcritical until a jump up to supercritical 
operation occurs; at this point, because of much lower bearing deflection, the rotor is well 
into the supercritical region. During run down, a larger phase angle change occurs before 
the jump than during run up, so that the phase angle just before the jump down is nearly 
equal to the 35 MN/m stiffness case. Above the critical speed, phase angle decreases for 
all stiffness cases, contrary to expectation. For the 35 MN/m stiffness case, phase angle 
begins to rise again with speed as the second critical speed is approached. 

The anomalous behavior of phase angle decreasing with increasing speed may be 
explained by looking at the difference in phase angle between the two bearings. This is 
shown in figure 12. At low speed, the phase angle difference is small, as expected for a 
translatory or bowed mode shape as shown in figure 3 (a) for the first critical speed. 

Above the critical speed, the two 
bearings are out of phase, as for 
the second critical speed of 
figure 3 (a). In the speed range of 
bistable operation (with variable 
bearing stiffness), phase angle 
during run up and run down 
shown in figures 11 and 12 may 
be studied along with the 
43 000 rpm mode shape plots of 
figure 9. During run up the bowed 
mode shape is nearly symmetric 
about the rotor axial midpoint; 

operation is still well within the 

subcritical region, resulting in small phase angles that are nearly the same for the two 
bearings. During run down, however, operation is supercritical at 43 000 rpm (above both 

the first and second critical speeds, as discussed above). The rotor has assumed a non- 

planar but near-conical mode shape, resulting in the two bearings being out of phase by 
90 degrees. Although not shown, the bearing at station 2 presents a more classical 



Figure 12. Phase angle difference between bearings 



Speed, rpm 

Figure 11. Station 10 phase angle 
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behavior, with phase angles increasing by nearly 180 degrees as the rotor passes through 
the critical speed. 

Effect of bearing damping. - It is known that ball bearings have very low inherent 
damping, but little data is available on actual damping coefficients. Lewis and Malanoski 
[10] mention values of 2.6 to 3.5 kN sec/m for a greased-packed bearing, but do not give 
the source of the data. Zeillinger and Kottritsch [11] describe experiments on an oil- 
lubricated 45 mm bearing. They took care to exclude damping of the bearing mount and 
include only the damping occurring within the bearing itself. Values of 0.5 to 1 kN sec/m 
were measured; for the smaller 25 mm bearing used in the present study one would 
expect less damping than this. Zeillinger and Kottritsch did note that additional damping 
usually occurs between the bearing outer race and its housing. However, this damping 
would be effective only for small vibration amplitude because of the small clearance 

For the data presented 
above, bearing damping of 
1.8 kN sec/m was assumed, 
figure 13 shows the effect of 
other damping values for the 
case of 12 g cm imbalance. 
Damping of 3.5 kN sec/m 
(twice that used above) reduces 
the peak amplitude considera- 
bly at bearing station 10 and 
separates the first and second 
critical speeds. There is still a 
region of bistable operation, 
however. A somewhat higher 
damping of 5.3 kN sec/m 
eliminates the bistable region 
entirely. This damping value is 
not large compared to what is prescribed for optimally damped rotors; however, it is 
more than can be expected from a ball bearing. Thus, as noted above, ball bearing 
supported rotors that are designed for operation near or above critical speeds usually 
incorporate dampers. 


between bearing and housing in typical installations. 



Speed, rpm 

Figure 13. Station 10 amplitude for various 
damping values; imbalance = 11 g cm 


CONCLUDING REMARKS 

Unbalance response data was presented for a rotor supported on ball bearings with 
accurate bearing stiffness calculated as a function of speed and load. Bearing stiffness 
was found to be a strong function of bearing deflection, with higher deflection producing 
markedly higher stiffness. Rotordynamic analysis showed that unbalance response varied 
nonlinearly with the amount of rotor imbalance. Moreover, the increase in bearing stiff- 
ness as critical speeds were approached caused a large increase in rotor and bearing 
vibration amplitude over part of the speed range compared to the case of constant bearing 
stiffness. Regions of bistable operation were possible, in which the amplitude at a given 
speed was much larger during rotor acceleration than during deceleration. A moderate 
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amount of damping eliminated the bistable region, but this damping is not inherent in ball 
bearings. 
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